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Abstract The dynamical analysis for coupled dark energy 
with dark matter is presented, where a complex scalar field 
is taken into account and it is considered in the presence of a 
barothropic fluid. We consider three dark energy candidates: 
quintessence, phantom and tachyon. The critical points are 
found and their stabilities analyzed, leading to the three cos¬ 
mological eras (radiation, matter and dark energy), for a 
generic potential. The results presented here enlarge the pre¬ 
vious analyses found in the literature. 

1 Introduction 

Observations of Type lA Supernova indicates that the uni¬ 
verse undergoes an accelerated expansion [1, 2], which is 
dominant at present times (^ 68%) [3]. Despite of ordinary 
matter, the remaining 27% is an unknown form of matter 
that interacts in principle only gravitationally, known as dark 
matter. The nature of the dark sector is still mysterious and 
it is one of the biggest challenges in the modern cosmol¬ 
ogy. The simplest dark energy candidate is the cosmological 
constant, whose equation of state wa = pa/Pa = is in 
agreement with the Planck results [3]. This attempt, how¬ 
ever, suffers from the so-called cosmological constant prob¬ 
lem, a huge discrepancy of 120 orders of magnitude between 
the theoretical prediction and the observed data. 

Among a wide range of alternatives, a scalar field is a 
viable candidate to be used. Its usage includes the canonical 
scalar field, called “quintessence” [4—8], and the scalar field 
with the opposite-sign in the kinetic term, known as “phan¬ 
tom” [9, 10]. Beyond the real scalar field case, a complex 
quintessence was also used in ref [11] to account the accel¬ 
eration of the universe. The t/( 1) symmetry associated with 
this complex scalar leads to a more sofisticated structure for 
the dark sector, and unless the stardard model of particle 
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physics is a very special case in the nature, there is no rea¬ 
son (apart from simplicity) not to consider a richer physics 
of the dark sector*. 

Another possibility of noncanonical scalar field is the 
tachyon and it comes from string theory. In the bosonic string, 
its ground state is the tachyon field, whereas in supersym¬ 
metric string theory a real tachyon is present in non-BPS 
Dp-branes, while the complex tachyon appears in a brane- 
anti-brane system [13]. The tachyon potential has a mini¬ 
mum [14, 15] and at this minimum the tachyon field behaves 
like a pressureless gas [16]. As soon as tachyon condensa¬ 
tion in string theory had been proposed, tachyon was also 
regarded as a dark energy candidate [17-19]. 

Still regarding the dynamical dark energy, there exists 
the possibility of interaction between dark energy and dark 
matter [20, 21], since their densities are comparable and, de¬ 
pending on the coupling used, the interaction can alleviate 
the coincidence problem [22, 23]. This approach was ap¬ 
plied to phantom and tachyon as well, in refs. [24-29]. 

When a scalar field is in the presence of a barothropic 
fluid (with equation of state Wm = Pm/Pm) the relevant evo¬ 
lution equations can be converted into an autonomous sys¬ 
tem and the asymptotic states of the cosmological models 
can be analyzed. Such approach was done for uncoupled 
dark energy (quintessence, tachyon held and phantom held 
for instance [30-34]) and coupled dark energy [21, 28, 35- 
41]. Since the complex scalar held and the coupled dark 
energy are generalizations of the real held and the unclou- 
pled case, respectively, we aim to study both possibilities 
together, in the light of the linear dynamical systems the¬ 
ory. Thus, we investigate in this paper the critical points that 
come from the evolution equations for the complex scalar 
held (quintessence, phantom and tachyon), considering the 


' A current example of a vector field that perhaps interacts with dark 
matter is the so-called “dark photon” (see [12] for a quick review). 
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possibility of interaction between the two components of the 
dark sector. The dynamical equations are derived and the 
critical point are analyzed, to find out what kind of universe 
can come up. As we have said, this is a natural extension 
of the previous works [21, 28, 37, 39]. As a result, we have 
found no new fixed points, however there are some crucial 
differences if compared with the case of real field. 

The rest of the paper is organized as follows. In section 2 
we present the basics of the interacting dark energy and the 
dynamical analysis theory, where we show the couplings we 
have used. In sections 3 and 4 we use the dynamical sys¬ 
tem theory for the canonical (and phantom) and the tachyon 
field, respectively. Within the respective section we show the 
critical points with their stabilities, and also the viable se¬ 
quence of cosmological eras (radiation-matter-dark energy). 
Section 5 is reserved for conclusions. We use Planck units 
{h = c = \ = Mpi = 1) throughout the text. 


2 Interacting dark energy and dynamical analysis 

We consider that dark energy is described by a scalar field 
with energy density and pressure p^, and with an equa¬ 
tion of state given by Wij, = p,j,/Pcj,. We assume that the scalar 
field is coupled with dark matter, in such a way that total 
energy-momentum is still conserved. In the flat Friedmann- 
Robertson-Walker background with a scale factor a, the con¬ 
tinuity equations for both components and for radiation are 

+3H{pij) + pij,) = ( 1 ) 

Pm + 3Hpm = (2) 

p,+4Hpr = Q, (3) 

respectively, where H = a/a is the Hubble rate, ^ is the cou¬ 
pling and the dot is a derivative with respect to the cosmic 
time t. The indices m and r stand for matter and radiation, 
respectively. The case of ^ > 0 corresponds to dark energy 
transformation into dark matter, while ^ < 0 is the transfor¬ 
mation in the opposite direction. In principle, the coupling 
can depend on several variables ^ = ^{p„t,p^,^,H,t,. ..), 
so we assume for the canonical scalar (quintessence) and for 
the phantom field the coupling is ^ = QPm^ [20, 21], while 
for the tachyon field the coupling is ^ = Qpmp^^/H [39], 
where 2 is a positive constant. The case with negative Q is 
similar and we will not consider it here because the minus 
sign of the case ^ < 0 can be absorved into the 0, instead 
of considering Q<Q. 

To deal with the dynamics of the system, we will de¬ 
fine dimensionless variables. The new variables are going to 
characterize a system of differential equations in the form 

2f' = /[2f], (4) 

where X is a column vector of dimensionless variables and 
the prime is the derivative with respect to log a, where we set 


the present scale factor ao to be one. The critical points Xc 
are those ones that satify X' = 0. In order to study stability of 
the fixed points, we consider linear perturbations U around 
them, thus X =Xc + U. Ai the critical point the perturbations 
U satisfy the following equation 

U' = /U, (5) 

where ^ is the Jacobian matrix. The eigenvalues of ^ de¬ 
termine if the critical points are stable (if all eigenvalues are 
negative), unstable (if all eigenvalues are positive) or saddle 
points (if at least one eigenvalue is positive and the others 
are negative, or vice-versa). 


3 Quintessence and phantom dynamics 


The complex scalar field <P can be written as ^> = 
where (j) is the absolute value of the scalar field and 9 is 
a phase. Both canonical and phantom fields are described by 
the Lagrangian 

^ = -V^[^d>^^dp0 + V{\<P\)), ( 6 ) 

where y(|^>|) is the potential for the complex scalar and 
we consider it depends only on the absolute value of the 
scalar field (j) = |^>|. We have e = 4-1 for the canonical 
field (quintessence) and e = — 1 for the phantom field. For 
a homogeneous field 0 = 0(f) and 0 = 0(f), in an expand¬ 
ing universe with Friedmann-Robertson-Walker metric with 
scale factor a = a{t), the equations of motion are 

e^ + 3eH(j)+V'{(j))-s(j)9^ = 0, (7) 

ee+(^3H+^^e=Q. (8) 

where the prime denotes derivative with respect to 0. In the 
uncloupled case Eq. (8) gives rise to the effective potential 
in (7), which is [11], where O) is an in¬ 

tegration constant interpreted as angular velocity. The first 
term in the brackets drives 0 away from zero and the factor 
may make the term decreases very fast, provided that 0 
does not decrease faster than . 

We assume the interaction between the scalar field with 
dark matter through the coupling Qpm^ and it enters in the 
right-hand side of Eq. (7). 

In the presence of matter and radiation, the Eriedmann 
equations for the canonical (phantom) field are 

= 0^ + y{<!>)+Pm+ Pr) , (9) 

77 = -^ (e0^ + e0^0^ + Pm + ^Pr) , (10) 


and the equation of state becomes 
_ P^ 0^-f0^0^-2£y(0) 

P0 02-4 0202 + 2ey(0)' ^ ^ 

We are now ready to proceed the dynamical analysis of 
the system. 
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3.1 Autonomous system 

The dimensionless variables are defined as 


XI = 




X2 = 

sfPr 


_ (i>e 


^ y=VM 


_ — yHr 5 = _i_ 'P — 

V3H’ V ’ - v/2 • 


( 12 ) 


The dark energy density parameter is written in terms of 
these new variables as 
Pi? 


^<l> = + exl +y^, 

so that Eq. (9) can be written as 

I I — 1. 


(13) 


(14) 


where the matter and radiation density parameter are defined 
by Qi — Pi/{3H^), with i = m,r. FromEqs. (13) and (14) we 
have that xi , X 2 and y are restricted in the phase plane by the 
relation 

0 ^ £x^ -f £x2 “fy^lj (15) 

due to 0 < £2(1, < 1. Notice that if y = 0 the restriction (15) 
forbids the possibility of phantom field (e = — 1) because for 
this case £2,i, < 0. 

The equation of state W||, becomes 




ev j + exl — 
exj+ex^+y^’ 


(16) 


which is a trivial extension of the real scalar field case. The 
total effective equation of state is 


^eff = 


I +Pr 


P(j> + Pm + Pr 


7 7 7 Z 

= £Xi + £x 2 — y ; 


(17) 


with an accelerated expansion for Weff < — 1 /3. The dynam¬ 
ical system for the variables xi, X 2 , x^, y, z and A are 


= —3xi +X 2 X 2 + 3^ey^A — ^eQ{l —xl — xl 



(18) 

dX2 , 

- =—3X7—X\X7X7—X7ri -, 

dN 2 1/72 

(19) 

^ Ls 

II 

1 

Who 

(20) 

dy v/6 

(21) 

dz_ u-idH 

dN ^ ^ dN’ 

(22) 


(23) 

where 


1 dH 3, 2 2 i\ 

(24) 


Point 

fii 

42 

43 

(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

(g) 

(h) 

9 

Q^-l 

-1 

“1 + ^ 

3x^ + \/bQx\ 

MQ 

2(A+Q) 

—3 + £A(A + (2) 

_3 

-1 

-1 

see the main text 
3(xf- 1) 

3(l+2fi) 

2 ( 2 .+e) 

0 

0 

0 

0 

Txi \J ^ ~ 

0 

0 
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Stability 

(a) 

3 

2 

1 

2 

saddle 

(b) 

Q(Q + 2.)+ 2 

e-+5 

unstable or saddle 

(c) 

2 

1 

saddle 

(d) 

2+ — 

1-^ 

saddle 

(e) 


see the main text 


(f) 

-5 \/6ai A 

D 2 

1 

unstable or saddle 

(g) 

P4J 

PSd 

saddle or stable 

(b) 

-3-h^ 

-2+^-y 

saddle or stable 

Table 2 

Eigenvalues and 

stability of the 

fixed points for the 


quintessence (phantom) field. 


3.2 Critical points 

The fixed points of the system are obtained by setting dxi /dN = 
0, dx 2 /dN = 0, dx-i/dN = 0, dy/dN = 0, dz/dN and dX /dN = 
0 in Eq. (18)-(23). When F = 1, A is constant the potential 
is y(0) = [30, 31]^. The fixed points are shown in 

Table 1 . Notice that X 3 and y cannot be negative and recall 
that £lr = z?'. Some of the fixed points do not exist for the 
phantom field because for those cases is negative. 

The eingenvalues of the Jacobian matrix were found for 
each fixed point in Table 1 . The results are shown in Table 
2 . 

The eigenvalues P 4 e and are 

_ 3(A+2e) 

P4e,5e — 4 {X+Q) ^ 

(25) 

At the first sight, one might think that the linear analysis 
would not give a complete description of the stability, be¬ 
cause all fixed points, but (f), have at least one eigenvalue 
equals zero. However, as pointed out in [42], fixed points 
that have at least one positive and one negative eigenvalue 
are always unstable, and methods such as center manifold 
[42] should be used to analyze the stability of the critical 
points that can be stable [((g) and (h)]. Even so, for almost 
all fixed points, but (a) and (c), X 3 = 0 , which means 0 — 00 . 

^The equation for A is also equal zero when = 0 or A = 0, so that A 
should not necessarily be constant, for the fixed points with this value 
of Xi. However, for the case of dynamical A, the correspondent eigen¬ 
value is equal zero, indicating that the fixed points is not hyperbolic. 
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Point 

Existence 

Xl 

X2 

X3 

y 

z 



12 

Weff 

(a) 

e = o 

0 

0 

any 

0 

0 


— 

0 

0 

(b) 

e = - 1-1 

-V6Q 

3 

0 

0 

0 

0 


1 

202 

3 

2Q^ 

3 

(c) 

any 

0 

0 

any 

0 

1 


- 

0 

1 

3 

(d) 

6 = -|-l 

-1 

V6Q 

0 

0 

0 


1 

W 

1 

1 

6^2 

1 

3 

(e) 

e = - 1-1 

2y/6 

3A 

0 

0 

2 V 3 

3A 


4 

A2 

1 

3 

4 

A2 

1 

3 

(f) 

e = -l-l 

any 


0 

0 

0 


1 

1 

1 

(g) 

any 

^/6 

0 

0 

/2e(fi+A)+3e 

0 


-e(e+2.) 

2(G+A)+3e 

-Q 

2{l+Q) 

V 2(A+(2)2 


Q(Q+2.)+3s 

(7+02 

X-\-Q 

(h) 

any 

eA 

V6 

0 

0 


0 


-1-f^ 

1 

1 

- 1 - 


Table 1 Critical points (xi, X2, X3, y and z) of the Eq. ( 18 ) for the quintessence and phantom field. It is shown the condition of existence, if any, 
of the fixed point (point (a) exists only for Q = 0 , for instance, while the point (b) does not exist for the phantom field.). The table shows the 
correspondent equation of state for the dark energy ( 16 ), the effective equation of state ( 17 ) and the density parameter for dark energy ( 13 ). 


However, this limit implies that X2 0 0 /H —00 as well, 
provided that H is finite. This issue occured for the points 
(b) and (d) to (h), as can be seen in the Table 1 whose math¬ 
ematical inconsistency indicates that the these critical points 
are not physically acceptable. 

Since all critical points are similar to those ones found in 
the literature [21, 28, 30] we reproduce the main results in 
the appendix 6 for the sake of completeness, which are valid 
for the case of real scalar field. 


3.3 Summary 

From the eight hxed points presented in the quintessence 
(phantom) case, only (a) and (c) are physically viable and 
they describe the sequence: radiation —matter. Both of them 
are unstable, however there does not exist a point that de¬ 
scribe the dark-energy-dominated universe. Thus, the extra 
degree of freedom due to the phase 9 spoils the physically 
acceptable hxed points that exist for the case of real scalar 
held, indicating that the dynamical system theory is not a 
good tool when one tries to analyze the complex quintessence 
(phantom). 


4 Tachyon dynamics 

The complex tachyon held <P — where ^ is the abso¬ 
lute value of the tachyon held and 0 is a phase, is described 
by the Born-lnfeld Lagrangian 

^Bi = -V^V{\cP\)^l-d^‘<Pd^<P, (26) 

where V(^) is the tachyon potential, which depends only on 
the absolute value of the scalar held 0 = |^>|. For a homoge¬ 
neous held 0 = 0(f) and 0 = 0(f), in an expanding universe 


with Friedmann-Robertson-Walker metric, the Lagrangian 
becomes 

ifB/ = -aV(0)^l-02_0202^ (27) 

where a = a(f) is the scale factor. The equations of motion 
for 0 and 0 are respectively 


1 — 0^02 — 02 

0 

1 _ 0202_02 


(1-1)292) 


3/20-1)92 _ 

1-<)292 + V{0) -^2, 


-3He + 


200 


0(1 — 0202 — 02 ) 


= 0 , 


(28) 

(29) 


where the prime denotes time derivative with respect to 0. 
When the phase 0 is zero, we recover the well-known equa¬ 
tion of motion for the tachyon held. 

The Friedmann equations for the complex tachyon, in 
the presence of barothorpic Huids, are 


H^ = - 


F(0) 


3 y-y/l — 0202_02 


+ Pm + Pr 


(30) 



/ V( 0 )( 0 ^-E 0 ^ 0 ^) 
y \/l — (l>^9^ — 0^ 


+ Pm + l^Pr 


The equation of state for the tachyon held yields 


W^ = ^ = 02 + 0202_l^ 

P(j> 


(31) 


(32) 


thus, the tachyon behavior is between the cosmological con¬ 
stant one (W|) = — 1) and matter one (w^ = 0). 

The interaction between the tachyon and the dark matter 
is driven by the coupling ^ = Qp^p^ip/H, which in turn we 
consider that it modihes the right-hand side of Eq. (28). With 
this form of coupling, the time dependence of the coupling 
is implicit in the Hubble parameter H. We are now ready to 
proceed the dynamical analysis of the system. 
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4.1 Autonomous system 


4.2 Critical points 


The dimensionless variables for the case of tachyon field are 


_ 1 




r = vv" 

-t — -pTT 


(33) 


■*1 = 0 , *2 = 00 , *3 = / 7 ^, 

^ 1 =_ YL 

y372’ 

Since 0 and 9 are dimensionless variables, (j) has dimension 
of time. 

The dark energy density parameter is written in terms of 
these new variables as 
„2 


£2(1) — 


P±_ 

3//2 






so that Eq. (30) can be written as 

£2(1) + £2,„ + Qr = 1 , 


(34) 


(35) 


where the matter and radiation density parameter are defined 
by £2i = Pi/{3H^), with i — m, r. From Eqs. (34) and (35) we 
have that xi , X 2 and y are restricted in the phase plane by the 
relation 


0 < *1 +*2 +y^ < 1 , 


(36) 


due to 0 < £2^ < 1. In terms of these new variables the equa¬ 
tion of state is 

w,j, =Xi+X 2 -l, (37) 

which is clearly a trivial extension for the complex scalar 
field. The total effective equation of state is 


^eff = 


P(l>+Pr 


= -y 


l-x]-xl + j, 


(38) 


“f Pm “t" Pr 

with an accelerated expansion for Weff < — 1 /3. The dynam¬ 
ical system for the variables x\, X 2 , * 3 , y, z and X are 


dx\ 

In 


dX2 

In 


= -(l-*l-* 2 ) X 

3xi-v^yA-f3e( l-z2- 




= -XiX2*3 -32C2(1 -Xj-X 2 ), 


dX 2 2 , *3 

— --X 1 X 3 + - 


2 I .2 3y2(l-;c2_^2) 

3 + Z- , =— 

1 2 2 

1 -X\-X2 


-|-X2X3,(39) 

(40) 

, (41) 


dN 2 


^ 2 3y2(l—X?—x,) 

-73xiyA + 3+z^- _ 4- I L 


1 - x^ - x^ 


dz ^ z 


2 3 y 2 (l-^ 2 _^ 2 ) 

3 H-Z — 


1 - Xj - X2 




(42) 

(43) 

(44) 


The fixed points of the system are obtained by setting dx\/dN = 
0, dx 2 /dN = 0, dx-i/dN = 0, dy/dN = 0, dz/dN and dX /dN = 
0 in Eq. (39). When F = 3/2, A is constant the potential has 
the form found in refs. [32, 43] {V{(j)) (j>^^), known in the 

literature for both coupled [28, 29] and uncoupled [17, 18] 
dark energy^. The fixed points are shown in Table 3. Notice 
that X 3 and y cannot be negative and recall that £2^ = z^- 
The fixed points yc, Xf and yy are shown below 


yc = 


WX^ + 36-X^ 


Q , \/e2+4 


(45) 

(46) 



The eingenvalues of the Jacobian matrix were found for 
each fixed point in Table 3. The results are shown in Table 
4. 

The points (al)-(a4) correspond to a matter-dominated 
solution, since £2,„ = 1 and Wgff = 0. They are saddle points 
because at least one eigenvalue has an opposite sign. The 
point (a4) is actually the point (al), with Q=\. Points (b), 

(c) and (d) are radiation-dominated solutions, with £2r = I 
and Weff =1/3. The difference between them is that (b) 
and (c) have Wij, = 0, while (d) has w,;, = — 1 and admits any 
value for X3. They are unstable [(b)] or saddle points [(c) or 

(d) ]. 

The point (e) is in principle a dark-energy-dominated so¬ 
lution with Q/j) = 1 and Wgff = w^ = — 1, whose existence 
is restrict to A = 0. However a careful analysis shows that 
the Jacobian matrix for this critical point has zero eigenvec¬ 
tor, thus it cannot be considered. Points (fl) and (f2) are also 
a dark-energy-dominated solution {Q,j, = 1) whose equation 
of state depends on A, which in turn can be either constant or 
zero. The case with constant A are shown in the Table 3 and 
an accelerated expansion occurs for A2 < 2/V^. For A = 0 
we recover the point (e). The eigenvalues /r 2 , P 4 and ps of 

the fixed point (fl) and (f2) are always negative. For these 

2^2 2 

points < 1, then the first eigenvalue is also negative if 
2 = 0, A < 0 or V3QXyc < 3. Therefore, the point (f2) de¬ 
scribes a dark-energy-dominated universe and can lead to a 
late-time accelerated universe if the requirement /Ti < 0 is 

^The equation for A is also equal zero when xi = 0, y = 0 or A = 
0 , so that A should not necessarily be constant, for the fixed points 
with these values of xi or y. However, for the case of dynamical A, the 
correspondent eigenvalue is equal zero, indicating that the fixed points 
is not hyperbolic. 
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Point 

Xl 

X2 

X3 

y 

Z 

W(l> 

^(j) 

Weff 

(al) 

±1 

0 

0 

0 

0 

0 

0 

0 

(a2) 

0 

±1 

0 

0 

0 

0 

0 

0 

(a 3 ) 

1 

0 

3 

2 

0 

0 

0 

0 

0 

(a 4 ) 

-Q 


0 

0 

0 

0 

0 

0 

(b) 

any 

±3/1-xf 

0 

0 

±1 

0 

0 

1 

3 

(c) 

1 

0 

2 

0 

±1 

0 

0 

1 

3 

(d) 

0 

0 

any 

0 

±1 

-1 

0 

1 

3 

(e) 

0 

0 

any 

1 

0 

-1 

1 

-1 

(fl) 

Aye 

0 

0 

yc 

0 

A^v^ 1 

3 ^ 

1 


(f2) 

AVe 

0 

V3lyc 

2 

yc 

0 

A^? 1 

3 ^ 

1 


(g) 

-Q 

0 

0 

0 

0 

G"-i 

0 

0 

(hi) 

Xf 

0 

0 

yr 

0 

Xy — 1 


x/y/X . 

Vi 

(h2) 

Xf 

0 

y/32.yf 

2 

yf 

0 



x/y/X . 

Vi 


Table 3 Critical points (xi, X2, X3, y and z) of the Eq. (39), for the tachyon field. The table shows the correspondent equation of state for the dark 
energy (37), the effective equation of state (38) and the density parameter for dark energy (34). 


Point fii 




/i4 Stability 


(al) 

6(1±G) 

(a2) 

0 

(a3) 

6(l-f0) 

(a4) 

0 

(b) 

6xf 

(c) 

6 

(d) 

-3 

(e) 

-3 

(fl) 

V^GAv,- 3 (i 

(f2) 

v^GAv,- 3 (i 

(g) 

-3(l-G") 

(hi) 

1 

1 

(h2) 



U 

6 

_3 

S 

j 

! 

? 

6(1-G^) 

3 ^ 

2 

1 

2 

! 

2 

6xf 

2 

2 

1 

-2 

-2 

2 

1 

-3 

2 

4 

1 

-3 

0 

-3 

-2 

-3(1-^) 

A^v^ 

2 

l^yl , 

2 

1 


A2y2 

X^yl , 


2 3 

2 

2 3 

2 2 

-3(1-G") 

3 

2 

3 

2 

1 

2 

-3(1-xy) 

3 

2 

n vT 

3 f 

2 1 ^3 

-^Axyyy-3(1 -xy) 

— ^Axyyy 

3 f x/yfX 
n vT 

3 / x/y/X 
2 V v/3 


saddle 

saddle 

saddle 

saddle 

unstable 

saddle 

saddle 

stable 

saddle 

stable for A < 0 or (2 = 0 

saddle 

saddle 

stable 


Table 4 Eigenvalues and stability of the fixed points for the tachyon field. 


statisfied. On the other hand, (fl) is a saddle point. The ef¬ 
fective equation of state depends only on A, so the coupling 
Q only changes the property of the fixed point. 

The point (g) is also a saddle point with a matter-dominated 
solution, however, different from (al)-(a4), the equation of 
state for the dark energy is no longer zero, but depends 
on Q, leading to an universe with accelerated expansion for 

< 2/3. For this point the coupling is restrict to values 
0<Q^<1. 

The last hxed points (hi) and (h2) are valid fox Xf 7 ^ 0^, 
for 2 7 ^ 0 and for constant A, and its behavior depends on 
Q. In order to have < 1, we must have Q > 0 for the case 
with plus sign in x/ (46), while we have Q < 0 for the minus 
sign case. We restrict our attention for the plus sign case. 
When (2 —^ 00 , X/ —> 0 and y/ —> 1, in agreement with the 

^The case for xy = 0 is the fixed point (e). 


restriction (36). In addtition, as pointed out in ref. [39], the 
fixed points exists for some values of A > 0 and Q, due to 
Eq. (36). Both fixed points have similar behaviour, however, 
(hi) is a saddle point, while (h2) is stable. Such difference 
is due to the eigenvalue (Table 4). The eigenvalues jX 4 
and jx$ are always negative because Wgff is between zero 
and minus one. The first eigenvalue is also negative because 
Q,j, < 1, thus Xy — 1 < — 1, therefore > -Vy, since 

Xf is always positive. Therefore, the point (h2) can lead to a 
late-time accelerated universe, depending on the value of A 
and Q. 

As in the case of quintessence and phantom, the hxed 
points that have X 3 = 0 [(al), (a2), (a4), (b), (fl), (g) and 
(hi)] indicate that 0 —3 00 and therefore X 2 = (j)9 —i' °° as 
well. However this limit is in contradiction to what is pre- 
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sented in Table 3 for xi showing that these seven critical 
points are not physically acceptable. 

All fixed points reproduce the previous results in the 
literature [28, 32, 39, 43] and they are generalizations of 
those analyses, with same stability behaviour for the criti¬ 
cal points. This indicates that the degree of freedom due to 
the complex scalar has no effect on the stability and on the 
evolution of the system of equations, when compared with 
the case of real scalar field. 

4.3 Summary 

The critical points showed in the tachyonic case describe 
the three phases of the universe: the radiation-dominated 
era, the matter-dominated era, and the present dark-energy- 
dominated universe. The matter-dominated universe can be 
described by the saddle point (a3). There are two points that 
can represent the radiation-dominated era: (c) and (d). The 
two points are saddle, with the additional difference that the 
point (d) has an equation of state for dark energy equals to 
minus one. 

A tachyonic-dominated universe is described by the point 
point (f2) and (h2). The point (f2) can be stable only if the 
coupling is zero or A < 0. The last fixed point (h2) is stable 
and can describe an accelerated universe depending on the 
value of A and Q. 

From all the critical points, the cosmological transition 
radiation —^ matter —dark energy is achieved considering 
the following sequence of fixed points: (c) or (d) —^ (a3) —^ 
(f2) [A dependent] or (h2) [Q and A dependent]. Although 
the sequence is viable, the form of the potential dictates 
whether the hxed points are allowed or not. Among several 
possibilities in the literature, the potential V((p) 0^", for 

instance, leads to a dynamically changing A (either if A —0 
for 0<n<2, orA—for n > 2) [19]. A dynamically 
changing A is allowed for the fixed points (a3), (c) and (d). 
On the other hand, points (f2) and (h2) require a constant A, 
implying F(([)) [17, 18, 28, 29]. 

5 Conclusions 

In this paper we studied coupled dark energy using a com¬ 
plex scalar field, in the light of the dynamical system theory. 
There were analyzed three possibilities: quintessence, phan¬ 
tom and tachyon field. All three possibilities are known in 
the literature for the real field [21, 28, 37, 39], and for un¬ 
coupled and complex quintessence field [33]. Thus, a natu¬ 
ral question that arises is how a complex scalar held changes 
the previous results and if there are new hxed points due to 
the complex held. Although some equations for the dimen¬ 
sionless variables are trivial extensions of the real held case 
(e.g. the equation of state for the scalar held), the differential 


equations were generalized. All hxed points found here are 
in agreement with the previous results, with no new hxed 
points, however there are some crucial differences. For the 
quintessence and the phantom there is a contradiction be¬ 
tween the hxed points X 2 and xt, when the latter is zero. This 
situation occurs for almost all hxed points and the only two 
exceptions are unstable points that represent respectively the 
radiation and matter era, so the dark-energy-dominated uni¬ 
verse is absent. Therefore the extra degree of freedom spoils 
the results known in the case of real scalar held. For the 
tachyon held all the critical points are also similar to the real 
held case, with same stabilities. Therefore, the extra degree 
of freedom due to the complex tachyon held plays no role 
on the stability of the critical points. Although the results 
presented here enlarge the previous results found in the lit¬ 
erature, with the generalization of the equations of motion, 
the dynamical system theory does not provide further infor¬ 
mation in what is already known for the case of real scalar 
held, letting open the possibility of studing complex scalar 
helds by other ways of analysis. 


6 Appendix 

In this appendix we reproduce the results in the literature 
[21, 28, 30] regarding the hxed points presented in Table 1, 
for the sake of completeness. 

The hxed point (a) is a saddle point which describes a 
matter-dominated universe, however it is valid only for Q = 
0. The other possibility of matter-dominated universe with 
arises from the hxed point (b). This point is called “0- 
matter-dominated epoch” (([)MDE) [21] and it can be either 
unstable or a saddle point. However, due to Q.^ = 2(2^/3 <C 

1 , the condition 2^ <C 1 should hold in order to the point be 
responsible for the matter era. Thus, fii and /i 2 negative, 
while /is is always positive and 114 is positive for Q{X + 
Q) > —3/2. Therefore (b) is a saddle point. 

The radiation-dominated universe is described by the crit¬ 
ical points (c), (d) and (e), only for the quintessence held. 
The hrst two points are saddle, as it is easily seen in table 

2 , and the last one had its stability described numerically 
in [21]. However, both (d) and (e) are not suitable to de¬ 
scribe the universe we live in, due to nucleosynthesis con¬ 
straints [44, 45]. The nucleosynthesis bound < 0.045 
[46] implies ^ > 3.7 for the point (d) and A^ > 88.9 for 
the point (e). Thus, the requirement for the point (d) is not 
consistent with the condition of point (b) and the constraint 
on A^ does not allow a scalar held attractor, as we will se 
soon. Therefore, the only viable cosmological critical point 
for the radiation era is (c). 

The point (f) is an unstable or saddle point and it does not 
describe an accelerated universe. The last possibility for the 
matter era is the point (g), with eingenvalues showed in table 



2 and Eq. (25). Since Wgff — 0 for |A | ^ \Q\, the fixed point 
is either stable or stable spiral, hence the universe would not 
exit from the matter dominance. 

On the other hand, the point (g) can lead to an accel¬ 
erated universe, for the quintessence field case (e = -fl), 
provided that 3 < A (A -f Q), because Q(j) < 1, and Q> Xjl, 
from Weff < —1/3. Regarding A > 0, the two eigenvalues 
/Ti and IJ .2 are always negative and since Q > 3jX — X, the 
behaviour of llAd, 5 d depends on the second term in the square 
root of (25) 

, ^ 8[3-A(A+0][3£ + 26(A+0] 

3 (A+22)2 ■ 

From the condition 3 < A (A + 2) we have A < 0, and if 
A < 1 the fixed point is stable. Otherwise, i.e. A > 1, the crit¬ 
ical point is a stable spiral. Thus, the value of the coupling 
dictates which behaviour the fixed point will have: stable for 
3/A—A <2<2*or stable spiral for Q> Q*, where 2* is 
the solution of A = 1. However, even in the case where one 
can get ~ 0.7 [45,47], there are no allowed region in the 
( 2 , A) plane corresponding to the transition from ())MDE to 
scaling attractor [21]. Thus, it is hard to gather the condi¬ 
tions for the point ())MDE and the point (g). For the case of 
the phantom field (e = —1), the condition / > 0 implies 
22(2 + ■^) >3. Hence, ^4 <0 and > 0, and (g) is a sad¬ 
dle point. 

The last fixed point (h) leads to an accelerated universe 
provided that A2 < 2. With this condition, the eingenvalues 
IJ. 2 , IJ .4 and Ids are always negative. The first eigenvalue /Xi is 
also always negative for the phantom field, and it is for the 
quintessence field with the condition A (A +2) <3. There¬ 
fore, the point is stable if the previous conditions are satis¬ 
fied. 
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